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On the Accuracy of the Numerical Solution of 
the Dirichlet Problem by Finite Differences ' 

J. L. Walsh 2 and David Young ^ 

This paper derives numerical bounds for the error, in certain closed regions, of the differ- 
ence analog of the Dirichlet problem. It is concerned only with the difference between the 
exact solution of the difference equation and the solution of the Dirichlet problem. The 
error bounds obtained involve quantities which can actually be computed, such as the mesh 
size, and the oscillation and modulus of continuity of the given function on the boundary. 
So far as the method is concerned, the chief novelty is the use of the difference analogs of 
harmonic measure and the Schwarz Alternating Process. 

1. Introduction 

Although finite difference methods afford a powerful tool for obtaining numerical solutions 
of partial differential equations, little is known about the accuracy. It is the purpose of this 
paper to derive numerical bounds for the error, in certain closed regions, of the difference 
analog of the Dirichlet problem. We shall be concerned only with the difference between the 
exact solution of the difference equation and the solution of the Dirichlet problem. The error 
bounds that we obtain involve quantities that can actually be computed such as the mesh 
size, and the oscillation and modulus of continuity of the given function on the boundary. So 
far as method is concerned, the chief novelty is the use of the difference analogs of harmonic 
measure and the Schwarz Alternating Process. 

Gerschgorin [4] ^ derived error bounds for boundary value problems associated with elliptic 
partial differential equations. These, and similar bounds derived by Collatz [2] and Mikeladze 
[10] involve bounds, in the closed region, of certain partial derivatives of the solution of the 
differential equation. However, the solution of the differential equation itself is not known, 
to say nothing of its derivatives, (although approximate values for the derivatives may some- 
times be found by examining the corresponding difference quotients). Also, it may happen 
that although the derivatives in question are not bounded in the closed region, the solution of 
the difference equation may still converge to the solution of the differential equation. 

Rosenbloom [13] presented an error bound for the Dirichlet problem which is closely 
related to Gerschgorin's but which utilizes special properties of harmonic functions. By use 
of well-known inequalities giving bounds for the partial derivatives of a harmonic function at 
an interior point in terms of the oscillation on the boundary and the distance from the bound- 
ary, he obtains upper bounds for the derivatives in closed subregions. Then, by solving the 
difference equation on subregions and, as the mesh size approaches zero, letting these sub- 
regions approach the given region, Rosenbloom obtains an error bound involving co*(5), the 
modulus of continuity in the closed region of the solution of the Dirichlet problem. However, 
Rosenbloom does not discuss the question of finding co*(5) in terms of co(5), the known modulus 
of continuity of the given function on the boundary. Furthermore, in practical numerical 
work one would not wish to change the boundary of the network as the mesh size is decreased. 

In section 3 we use the explicit solution of the difference analog of the Diriclilet problem 
for the rectangle, obtained by Le Roux [9] and by Phillips and Wiener [12], to derive an error 
bound involving bounds for the derivatives of the given function on the boundary. Wasow 
[15] has obtained error bounds for the rectangle which are applicable when the boundary values 
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have bounded third derivatives. Our results are somewhat more general in that we allow the 
second and third derivatives to fail to exist at a finite number of points. Also, we obtain error 
bounds for the case where only the second derivative is bounded, and we allow the first and 
second derivatives to fail to exist at a finite number of points. In sections 4 and 5 harmonic 
measure and its finite difference analog, discrete harmonic measure, are used to obtain bounds 
for the moduli of continuity of harmonic and discrete harmonic functions, respectively, in 
certain regions in terms of a)(5). Phillips and Wiener [12] showed the existence of such bounds, 
whereas we obtain suitable bounds in a precise numerical form. Then, in section 6 these upper 
bounds are used to yield a uniform error estimate for the closed rectangle in terms of the oscillia- 
tion on the boundary and 5(co). In section 7 the extension to regions made up of two or more 
overlapping rectangles is discussed. 

2. Discrete Harmonic Functions 

Let h and k be arbitrary positive numbers, and let L[h,k] denote the set of points {x, y) 
such that both xjh and yjk are integers. Let 12 be a simply connected closed region with interior 
R such that the boundary aS of 12 consists of straight lines, each of which is parallel to a co- 
ordinate axis and contains a point of L[h,k]. Let 12^, denote the subset of points of L[h,k] 
contained in 12. Two points (xi,yi) and (0:2,2/2) of L[h,k] are adjacent if 

[{xi-x,)/hY+[(yi-y2)/kY=i. (2.1) 

A point of Ql is an interior point of 12z, if the four adjacent points belong to 12^^. All other points 
of 122:, are boundary points. We let R^ and Sl denote, respectively, the set of interior and 
boundary points of ^l. Evidently, we have Rl^R and Sl^S. 

A function U(x,y) defined on Ql is said to be discrete harmonic,^ (d. h.) in Rl if it satisfies 
the difference equation 

S'[Uix,y)]^[2aVil + <.')][U(x+h,y) + U{x-h,y)-2U{x,y)] 

+ l2lil + .r')][U(x,y+k) + U{x,y-k)-2Uix,y)]=0, (2.2) 
where 

<T = k/h. (2.3) 

The finite difference analog of the Dirichlet problem is the following problem: given a 
function f(x,y) defined on S, to find a function U{x,y) defined on Qj,, d.h. in R^ and coinciding 
with f{x,y) on &. The existence and uniqueness of a solution of this problem for bounded 
regions is easy to prove; see for instance Gerschgorin [4]. The convergence to the solution 
of the Dirichlet problem has been proved using nonconstructive methods by Le Koux [9], 
Phillips and Wiener [12], and others. 

3. Error Estimate for the Rectangle Under Differentiability Assumptions 

Let 12 be bounded by the lines x=0, x=a, y=0, and y=b, where a=Ah, b=Bk and A and 
B are positive integers. 'Letf(x,y) be defined and continuous on S and let u(x,y) and U(x,y) 
denote respectively the solution of the Dirichlet problem and its finite difference analog with 
boundary values determined hyf(x,y). In this section we shall derive an upper bound for the 
error U(x,y)—u{x,y) in the region Q^, under certain assumptions about the derivatives oif{x,y). 
lifi(x,y) is defined and continuous on S, we denote generically by Ui(x,y) and Ui(x,y) the solu- 
tions respectively of the Dirichlet problem and its finite difference analog with boundary values 

First, one can verify directly that the function Ui(x,y) defined by 
«i (x,y) =/(0,0) + [/(a,0) -/(0,0)] (x/a) + [f(0,b) -/(0,0)] (y/b) 
+[f(a,b)+fiO,0)-J{a,0)~J{0,b)]ix/a){y/b), (3.1) 

» Heilbronn [5] introduced the term discrete harmonic function, and studied the properties of these functions. 
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which is hnear on every Hne parallel to either coordinate axis, is both harmonic in U and d.h. in 
Rl. Thus i{fi(x,y) is defined on S and equal to Ui{x,y) on S, then Ui(x,7j) = Ui(x,y) in Q^. 

Next, ii f2(^,y) =f(^,y) —fi(^,y) on S, then we have by linearity of Laplace's equation and 

(3.1) 

U(x,y) - u (x,y) = [ U2 (x,y) + Ui (x,y) ] - [u2 (x,y) + u^ (x,y) ] 

= U2{x,y)~U2(x,y). (3.2) 

Evidently both U2{x,y) and U2{x,y) vanish at all the corners of Q. 

We study arbitrary boundary values J{x,y) by studying in turn four functions each of 
which has the boundary values y2(^, 2/) ^^ ^^^ side of 12 and values identically zero on the other 
three sides. By the symmetry of the situation, it is sufficient to study in detail only one of the 
latter functions. Let us set J^(x,y)=f2(x,0) when y=0 and j's(^,y)=0 elsewhere on S. Since 
J2(x,y) is continuous and vanishes at the corners of U,f^(x,y) is continuous on S. 

It can be verified directly ^ that if the Fourier series of /a (a;, 0) converges tofz{x,0), then we 
have 

^3(x.y)=2 A„ sm in.x/a) ^^^^Uf,^/"' ' (3-3) 

where ^ 

An = (2/a) j /3(^0) sin {mrtla)dt, (3.4) 

A:=(2h/a)j: sin {nirjhla)J,{jhfi), (3.6) 

;=i 

and where m and n satisfy the relation 

sinh (m7rZ:/2a) = 0- sin {mrhl2a). (3.7) 

Now let g{x) =f^(x,0). We now prove the following theorem: 

Theorem 3.1 Let g(x) satisfy the following conditions: (a) g{x) and its first (s—2) deriva- 
tives are continuous ^ 0<x<a; (b) g^^~^\x) fails to exist or fails to be continuous for, at most, a 
finite number of points in the interval 0<x<a, and g^^~^^(x) is uniformly continuous in each open 
interval in which it exists and is continuous; (c) g^^\x) exists, except possibly for a finite number of 
points, and is bounded, 0<x<a. 

Then uniformly for all (x,y) eQ^, "^e have 

\U,(x,y)-u,(x,y)\<(2 + C)D2(h/a), (3.8) 

if s = 2, and 

1 U,(x,y) -u, {x,y) ! < (1 + Ce)D,{hla)\ (3.9) 

if s = Z, where 

(7=(l/24)7r2(l + (T2) coth {awb/a)/ae, (3.10) 

D2=2{a/Ty[2j,Mi/a+M2], D^=2{alTrY[2{J2+l)M2la+M,] (3.11) 

a:-(2/7ro-) sinh"^ (cr), (3.12) 



and that ' 



where 



• See, for instance [8], pages 95 and 96. 

7 See, for instance, Le Roux [9] or Phillips and Wiener [12]. 
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and Jij i—1^2, denotes the number oj points in the open interval 0<Cx<Ca at which g^^^(x) does not 
exist or is not continuous . Here Mi,i=l,2,3, denotes the least upper bound of the modulus ofg^^^(x). 

Proof. 

Lemma 3.1. Ifm and n satisfy (3.7)^ then 

0<7i-m<(l/24)(H-a-2)7rW(A/a)2, (n=l,2, . . .,A). (3.13) 

Proof. Let p=m7rk/2a, q=nTh/2a. Then (3.7) takes the form 

sinh ^ = 0- sin g. (3.14) 

We now study the function p(q) defined by (3.14) in the interval 0<g<7r/2. One can 
verify directly that p(0)=0, p^(0) = (r, and 

--^(l+O^m^ <o 



^ w)"~ n _L ^2 -^2^X5/2 



-(r(l+o-^)(l — 2o-^sin^g) cos q 
(l+cr^sin^g)^ 

where primes denote differentiation with respect to q. By the extended mean value theorem, 
since p^^(q) <0, we have p<aq. Therefore, m<n. On the other hand, we also have 

. //// ^ I <r <^(l + <^')|l-2cT^sin^g| 2 

'^ \^/i— (1 -\- a^ siir qY'^ 
since \(l-2a'sm\){l + a'sm\)-'^'\<l. 

Again using the extended mean value theorem, we get 

and 

\p(q)-aq\<{mqMl + cr'), 
and the lemma follows. 

Lemma 3.2. If m and n satisfy (3.7), and if a is given by (3.12), then 

m>an, {n=l,2, . . .,A). (3.15) 

Proof. Defining p and q by (3.14), we have 

j^(piq)lq) = [p'iq)q-p{Q)]lq'. 

By the mean value theorem, we have p{q)=qp\^) (0<^<g). Since, as shown in the proof 
of lemma 3.1, p"(g)<0, we have p'(q)^p'(^) and p^(q)q—p(q)<0. Therefore, the ratio 
p{q)lq is a nonincreasing function of g, (0<g<7r/2). Its minimum value in the interval 
0<q<'7r/2 is assumed when q=Tr/2. Therefore, p (q) / q'> p (it /2)/(t 12) = (2 /T)sijih.~^((T), and 
the lemma follows. 

Now for convenience, let us define the function 

smh[(jn.bla)il-y/b)l 
^^^ smh(m7r6/a) ^ 

We shall study its behavior as a function of both m and y, where m is assumed to be a con- 
tinuous variable. 

Lemma 3.3. If m> 0, then for <y<b we have 

> -v— Tm{y) > —{Try /a) coth (mirbja) exp [—miry /a], 
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Proof. Differentiating (3.16) with respect to m we obtain 

j^ Tm{y)=(7rb/a){(l — y/b) cosh [(m7r6/a)( 1 — ^/6)] sinh {rrnrhla) 

— cosh {mwb/a) sinh [(^7r6/a)(l— ?//6)] } sinh~Xm7r&/a) 
= (7rb/a)Tm{y)[(l—y/b) coth [(mwb/a){l — y/b)] — coth {rmrb/a)}. 
We note that for K>0, x>0, 

j^ {x coth {Kx))=coi}i {Kx)-Kx csch^ {Kx)=^ ^^" sinh^-g^)^ ^ * 

But since ^m\i{2Kx)'>2Kx, the last expression is nonnegative; therefore by the mean value 
theorem we have 

(l-?//6)coth[(m7r6/a)(l-?//6)]<coth(m7r6/a), (0<?/<6), 
and 

ir.(,)<o. 

On the other hand, 



dm 



= (,>Wcoth(.aw{(,-,/» -K2,Wa-^/>>' -r.(,)j 

= {irbla) coth (m7r6/a)S(l-2//J) . , , smh (m7ry/a) (y/6) r,(7/)l 

^ ' ^ ^ / /^v ;// ^ smh (m7r6/a) cosh (m7r6/a) ^^' ^ ^^^) 

For 0<y<b, the terms in the brackets are, respectively, nonnegative and nonpositive, whence 

d 



dm ^'«<^2/)>-(^2/M) coth {rmrbla) r^(t/). 



Finally, we note that 



provided^ <y<b. The lemma now follows. 
Lemma 3.4. If m and n satisfy (3.7), then 



\[-M^-^y^]m = m, 



< [coth (aw b I a) I ae]n-\ (1<^<^X (m,<mi<n,) (3.17) 



where a is given by (3.12). 

Proof. It is easily verified that if 0<<^, 0<K, then \te-^'\<(Ke)-K Therefore, 
y exp{-—m,7ry/a)<[{mir/a)e]~^. The lemma follows from lemmas 3.2 and 3.3, and from the fact 
that coth(a;) is a decreasing function of x^ for x>0. 

Lemma 3.5. If m and n satisfy (3.7), then 

\rm(y)-Tn(y)\<Cih/ayn', (0<y<b), {l<n<A), 
where C and a are given by (3.10) and (3.12), respectively. 

8 This inequality was proved by Phillips and Wiener [12j. 
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Proof. By the mean value theorem we have 

|r.(^)-r,(2/)i<Max^<^^<^|[^ rUv)]^^^^ 

The lemma follows from (3.13) and (3.17). 

Lemma 3.6. If An is defined by (3.4), then 

\An\<Ds/n\ (71 = 2,3), 



m— n 



(3.18) 



where Ds is defined by (3.11). 

Proof. The proof involves the use of repeated integration by parts of the integral ex- 
pression for An, and is similar to that given by Jackson[8, p. 13-14] for the case s = 2. We 
omit the details. 

We now define the function /4(a;,2/) : 

A 

f4{Xjy) = unless y=0 and fA(x,0) = y^,An sin irnrxla). 
For ^=2,3 and 0<x<a, we have 

|/4(X,0)-^(X)| = 



y^, ^n sin (nirx/a) 

n = A-\-l 



<Ds i] n-'<Ds( t-'dt=DsA'-'/(S'-l), 

n = A-\-l J A 



or 



\Uxfi)-g{x)\<[D,l{s-l)]{hlay-\ 



(3.19) 



On the other hand, if we replace fzixfi) by J^ixfi) in (3.4) and (3.6), we get A^=An 
(/i=l,2,...,^-l). Therefore, for {x,y) e L lh,k], we have, since sin (ATx/a) = Oy 



and 



Ui(x,y)—U4(x,y) = J2 ^n sin (mrx/a)[Tmiy)—Tn(y)]y 

n=l 



\U,{x,y)-u,(x,y)\<j:\An\-\r^{y)-TM\<CD,(h[ayj:n'- 

W = l 71 = 1 



by lemma 3.5 and (3.18). It can be shown that 

A = {hla)-\ {s=2) 

+\ogA = l-log{hla), (s = 3). 

Therefore, if 5=2, 

I U^ {^jV) — '^4 {x,y) I < C©2 {hi a) . 



A-l PA ( 

n=l Ji Q. 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



Evidently (3.22) could be used to obtain an error bound for the case s=3, and this was done 
in the original manuscript of the present paper. Since then, however, a paper by Wasow [15] 
has appeared which contains a stronger result for the case s==3 than the original form of theorem 
3.1. A combination of his methods and the original ones now yields an improvement (the 
present theorem 3.1) of Wasow's result in this case since we allow g^^^{x) and g^^^(x) to fail to 
exist for a finite number of points. 

Thus following [15], we obtain from lemmas 3.1, 3.2, 3.3, and the mean value theorem 

\Tm(y) — Tn(y)\<Cae(hlay(Tryla)n^ex]^(—anTyla). 
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By (3.20) and lemma 3.6 we have 

\U,{x,y)-u,{x,y)\<CD^ae{hla)\iryla)Y:^^V{-otniryla)<CeD^{hla)\ (3.24) 

00 

since for ^/^O^S^^p (— an7r?//a) = [exp {aTTyla)—\]~^<{aTryla)~^. 

71 = 1 

Since d.h. functions possess mean value properties, they assume their maximum and mini- 
mum values on the boundary as do harmonic functions. Since the functions [uz{x,y)—Ui{x^y)] 
and [Uz{x,y) — U^{x^y)] are harmonic and d.h., respectively, we have 

\u^{x,y)—u^{x,y)\<M2iX \g{x) — f^{x,())\, \Uz{x,y)—U^{x,y)\<M2ix \g{x)—f^{x,())\. 

0<x<a 0<x<a 

Therefore, in ^l 

\Uz{x,y)-u^{x,y)\<\Uz{x,y)-U^{x,y)\ + \U^{x,y)-u^{x,y)\ + \u^{x,y)-u^^^ 
<\U,{x,y)-u,{x,y)\ + 2 Max \g{x)~f,{xM- 

0<x<a 

By (3.19), (3.23), and (3.24), the theorem follows. 

4. Harmonic Measure 

Let 12 denote a simply connected region with interior li and boundary S. Let >S'' denote a 
subset of S consisting of a finite number of connected subsets of S. We define the harmonic 
measure, (h.m.), H[(x,y),S\Q] as the unique function which is harmonic and bounded in R, is 
continuous in 12 except perhaps at a finite number of points of S, and equals unity on S^ and 
zero on S—S\ The properties of harmonic measure have been studied in considerable detail, 
see Nevanlinna [11], chapter III. 

By analogy we define discrete harmonic measure, (d.h.m.), for regions of the type described 
in section 2, as follows: 2?z,[(a:,?/),iS",12j is a function d.h. and bounded in El, equal to unity on 
SLdS' and to zero on (S-S')nSL. 

For bounded regions the existence and uniqueness of h.m. is well known, see [11]. The 
existence and uniqueness for the half plane can be proved by the use of conformal mapping. 

The existence and uniqueness of d.h.m. for bounded regions follows from the existence 
and uniqueness of the solution of the difference analog of the Dirichlet problem. Later we shall 
prove existence and uniqueness of d.h.m. for a half plane and for certain other unbounded regions. 

We list some elementary properties of d.h.m., which are analogs of well-known properties 
of h.m. 

<HL[ix,y),S',QL] < 1 for all S'^S. (4.1) 

If S'^S^'^S, then 

Hd(x,y),S\Qr.] <Hd{x,y),S^\Qr.l (4.2) 

If aS"^ is included in the boundary of both Q^ and 12* where 12^^12* and if {x,y)eQ>L, then 

Hd{x,y),S\^j] <H^[{x,y), 8^,^11 (4.3) 

The first property follows at once from the fact that the maximum and minimum values 
of d.h. functions are assumed on the boundary. The second follows from this fact and from 
the fact that the expression [HL[{x,y),S'\9.i\—HL[{x,y)yS^j^L]} is nonnegative on Sl- The 
third follows since by the maximum and minimum principles HL[{x,y),S' j^l\>0 and hence 
{HL[{x,y),S',^l]—HiX{x,y),S',^L^>Q{oT any point {x,y) on the boundary of 12£,. This is the 
so-called principle of gehietserweiterung. 

We shall use h.m. and d.h.m. for two purposes. In this section and in section 5 lower 
bounds will be derived to enable us to obtain upper bounds for the modulus of continuity of 
harmonic and d.h. functions in a closed region in terms of their moduli of continuity on the 
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boundary. In section 7, upper bounds will be derived which will enable us to use the Schwarz 
Alternating Process and its difference analog for overlapping rectangles. 

Let U denote the rectangle of section 3, and let S^ denote the side contained in the line y=0. 

Theorem 4.1. // (x-a/2y+y^<d', then 

H[{x,y),S^,^]>l-vib/a)(8/a), (4.4) 

Hd{x,y),S',Qr.]>l-v{bla){8la), (4.5) 

where 

p{bla) = [4: + Tr'cot]i'(irb/a)yf'. (4.6) 

Proof. Let u{x,y), U(x,y) be harmonic and d.h. functions respectively in R vanishing 
on S—S^ and equal to sin(7rx/a) on S\ By (3.3) to (3.7) we have 

«(a:,y)=sin(.x/a) 5iHM?Ma^l^M, and ?7(^,y)=sin(.x/a) ^^^^[^"^"^/'^^^/'^/''^l 
^ ' ' sinh(7r6/a) ^ '^ v / / smh(m7r&/a) 

where sinh(m7rA:/2a)=cr sm('jrh/2a). 

^j lemma 3.1 we have m^l. Also by lemma 3.3 we have (<i/c/m)r^(?/)<0. Therefore, 
r^ (y) > Ti (y) , and we get 

U(x,y)>u{x,y). (4.7) 

Now, for 0<^<7r/2 we have ^ sin 6>(2/T)d', and hence for 0<x<a/2, we have 

sm{wx/a)>2x/a=l — (2/a) {a/2—x). But since sin [(7r/a)(a/2 + p)]=sin [{7rfa){a/2 — p)] for all p 
it follows that sin(7r^/a)> 1 — (2/a)|a/2— a;|, (Q<x<a), 



Also 

d 



^sinh[(7r6/a)(l-2//6)] 



= |(7r/a)cosh[(7r6/a)(l —y/b)]\ <(7r/a)cosh(7r6/a), 



and by the mean value theorem 

sinh[(7r6/a) (1 —y/b)] >sinh(7r6/a) — (7r?//a)cosh(7r6/(x). 
Therefore, 

1 -u{x,y) < 1 - { [1 - (2/a) \a/2-x\][l - (7r?//a)coth(7r&/a)] } < (2/a) \a/2-x\ + (Ty/a)coth{Trb/a). 

By the Schwarz inequality the last expression does not exceed 

{l/a)[(al2-xy-i-yrHb/a) < m)v{bla), 
provided (x-'a/2y+y^<5^. 

We now observe that H[{x,y) ,S\Q]> u(x,y) , HL[{:x,y),S' ,Qi\> U(x,y)^ The theorem now 
follows from (4.7). 

Now let ftx, denote the subset of L[h,k] such that y> and let El and Sl denote the interior 
and boundary of ^l, respectively. 

Theorem 4.2. For the region Ql, d.h.m. exists and is unique. 

Proof. Consider the sequence [Ht^{x,y)], where 

Ht'(x,y)=HL[(x,y),S^[]Sj;:\^t'l 

and where ^l^ denotes the rectangle with vertices (nh,0), (—nhfi), {nh,nk) and {—nh,nk). 
Evidently HL\x,y) exists, since d.h.m. exists for bounded regions. By (4.3) it follows that the 
sequence {HL^{x,y) ] is nondecreasing and bounded above by unity. Therefore, a unique limit 
exists which we denote by HL(x,y) . 

• This is Jordan's Inequality; see Copson [3, p. 136]. 
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Since HL\x,y) is d.h., it satisfies (2.2). Taking limits of both sides of (2.2), we find that 
HL(x,y) satisfies (2.2) and is therefore d.h. Evidently HL{x,y) <l, HL(^x,y) = l on S'CiSl and 
HL(x,y)=0 on S^ — Sl. Hence nL(x,y) has all the properties required for d.li.m. 

In order to establish uniqueness, we prove 

Lemma 4.1. // U(x,y) is d.h., bounded in Rl and vanishes on Sl, then U(x,y) vanishes in Rl. 

Proof. By hypothesis there exists a constant P such that for (x,y)€QL, we have | U(x,y) \ <P. 
Given any {x,y)eUL and any e>0, let a and b be chosen so that {a,b)eL[h,k] and 

a>P(x'+y'y^h{l/2)/2e, b>a. 

Let 12* denote the rectangle with vertices at (a,b), (—a,b), (a,0), and (— a,0), and let la denote 
the interval (— a<^x<^a). Since \U(x,y)\<P on aS* and vanishes on /„, we have 

\Uix,y)\<PH,S(x,y),S*~IayL]. 
By theorem 4.1 we have, replacing a by 2a, 

Hd(x,y),Ia,nn>l-vib/2a){8/2a), 
provided x'^-\-y^<d^. But by hypothesis, we have (x^+y^)^^^<2ea/Pv(]/2). Substituting, wo 

Hdix,y),Ia,nn>l-v{b/2a)e/Pv{l/2). 
Since b> a, we have v{b/2a) <v(l/2) ; hence, 

Hd{x,y),IaX]>^-elP, and \U(x,y)\<e. 

This proves the lemma. 

The uniqueness can now be proved by assuming two d.h.m.'s and showing that their 
difference vanishes identically. 

The proof of theorem 4.2 is complete. 

We note that the existence of a unique bounded solution of the difference analog of the 
Dirichlet problem for fiz, with bounded boundary values is almost immediate. If the boundary 
values are determined by g{x), then the limit of the absolutely and uniformly convergent series 

X) HL[(x,y) , fxh,UL]g(l^h) is d.h. in Rl, bounded in Ql and equals g(x) on Sl. Thus a solution 

exists. The uniqueness follows at once from lemma 4. 1 . 

We next consider the region x>0, y>0. To find d.li.m. for subsets of the line y=0 we 
perform a sign-changing reflection about the line x=0 and use theorem 4.2. Since the d.h.m. 
is zero on the line x=0, this can be done in such a way that the new function will be d.h. in 
Rl. Similarly to find d.h.m. for subsets of the line y=0 we reflect in the line x=0. 

By similar methods the existence and uniqueness of d.h.m. can be established for other 
regions such as the semi-infinite strip 0<x<a, y>0. 

5. Modulus of Continuity 

The modulus of continuity of a function /(x,?/) in a closed region ^ is defined by 

c.{8)=LUB\f{x„y,)-J(x2,y2)i 

where {xi,yi),{x2,y2) el2 and (o^i— ^2)^+ (2/1—2/2)^ <5^ Evidently if 5 is not less than the diameter 
of 12 then oo{8) equals the oscillation oif{x,y) in 12. 

For harmonic functions we prove the following theorem: 

Theorem 5.1. Let 12 denote a bounded simply connected closed region with interior R and 
whose boundary S is a closed Jordan curve with the following property: there exist constants 
ro^O and d>0 such that jor any point P of S there exists a circular sector with vertex at P, with 
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radius Tq and included angle 6 containing no point of R, Let u{x,y) he harmonic in R, continuous 
in U with modulus of continuity co(5) on S. If for some positive number D we have 

8<M^x{ro,D,D (ro/Dy-^^^^-} 

then the modulus of continuity of u{x,y) in 12 satisfies the inequality 

co*(6)<co[Z)(3/D)-/^-+^>] + (4M/7r)(5/Z))-/^-+^\ (5.1) 

where i^=Max (27r— ^,7r) and 

M=Max [u{x,y)]—Mm [u(x,y)] (5.2) 

(:c,y)eS (z,y)eS 

is the oscillation of u {x,y) on S. 

Proof. 

Lemma 5.1. Let Ir denote the interval —r<x<r, y=0 and let Cp denote the region 
x'+y'<p', 2/>0. // (x,y)eC,, then 

H[{x,y)Jr,y>0]>l-2p/irr. 

Proof. It is eas}^ to show that if <^ is the angle at (— r,0) from /;. to the circle C through 
the points (— r,0), (r,0) and (x,y), then H[(x,y),Ir,y'^0] = l — (t)/T. 

Now if p/r= (1 — cos </))/sin <^, then Cp is contained in the region bounded by C and the hne 
y=0. Hence, by the minimum principle for harmonic functions we have for {x,y)€Cp 

H[ (x,y) Jr,y > 0] > 1 - (2/x) tan-i (p/r) > 1 - 2p/7rr, 

and the lemma is proved. 

Lemma 5.2. Let Cr,,p denote the region x^-]-y^<r^, < tan"^ (?//;r) <\l/<27r and let Br,yp denote 
the bounding radii of Cr,^p. {If yl/=2Tr, then Br,rp denotes the line 0<x<r, y=0). If p<r and if 
{x,y)eCp,4,, then 

H[{x,y),Br. ,,a. ^] > 1 - (4M (p/ry*. 

Proof. We first consider the case \l/=7r. If x^-{-y^=r^ and y^O, then as in lemma 5.1 
we have II[{x, y), I r,y'^0] = 1/2. Hence, as we verify at once, 

H[(x,y),IrA] = 2H[(x,y)Jr.y>0]-h 

If {x,y)eCp,r, then by lemma 5.1 we have 

H[{x,y),Ir,Cr]>-^-W^)(p/r). 

The lemma can now be verified for the general case by mapping Cr,Tr onto Cr,^p by means of the 
conformal transformation {wlr) = {zlry'^. We omit the details. 

Now by the maximum principle for harmonic functions, for given 5 the maximum value of 
\U{xi,yi)—U{x2,y2)\ for (^^i— 0:2)^+ (y 1—2/2) ^< 5^ occurs when either (xi,2/i) or (^2,2/2) belongs to S. 
If {xi,yi)eS, then for all r such that 0<^r<ro there exists a circular sector 6^^ containing at 
least one point of R with vertex at (a?i,^i) with radius r>5 and angle xf/, such that Br,^p the union 
of the bounding radii is disjoint from R. Since the theorem is trivially true if {x2,y2)^S, we 
assume {x2,y2)^R' Let 12i denote the closure of the connected component of ^,2/2) for the region 
Qf]Cr,4.. Evidently ^icCr,^. Since for (x,y)e^if]S, we have 

H[(x,y),Br,^A.i]<^=II[(x,y).^r(]SM 

and for (x,y) contained in Qi and on the arc of C^,^, we have 

H[{x,y),Br,^,CrJ = 0, 
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then for all (x,y)e^i it follows that 

But by lemma 5.2 we have 

H[{x,,y2) A ^A, ^] > 1 - (4/7r) (6/ry^^ 

since the lemma is obviously true if we rotate the sector. 
Therefore, we have 

\u(x2,y2)-u(x,,y^)\<c,{r)H[{x2,y2)A.^f)SM]+MH[(x2,y2),S^ 

Now by the assumptions on 8 we hsive D{dlD)''^^'''^'^^ <D{ro/D)=rQ. Hence, we can choose 
r=Z>(5/Z))''^^^+^^ <ro, and the theorem follows. 

We now consider the case where 12 is an arbitrary simply connected region. We first prove 
the following general theorem, which is essentially a formulation for harmonic functions which 
is equivalent to Caratheodory 's form of the theorem of Milloux for analytic functions : 

Theokem 5.2. Let G be a Jordan subregion of |^Kl, whose boundary consists of a Jordan 
arc ai which passes through z=0 plus an arc a2 of T:\z\ = l. Let u{z) be harmonic and bounded in 
G, continuous in the corresponding closed region except at the end points of ai and a2, egual to unity 
in the interior points of ai and to zero in the interior points of 0:2. Then in every point z of G we have 

'i^(2)>l-(4/7r)tan-i|^|^/l (5.3) 

Proof. Let v{z) be conjugate to u{z) in G, and lai f{z)=QiX^[—u{z)—iv{z)]. Except 
perhaps at the end points of the arcs, on ai we have \f{z)\ = e~^^ and on 0:2 we have \f{z)\ = l] on 
these open arcs 1/(^)1 is continuous in the two-dimensional sense. Then by the form of Mil- 
loux^s theorem presented by Caratheodory, [1], 354, we have for zeG 

\o^\f{z)\ = -u{z) <-l + {4.lir)t^n-'\z\'i^ 
as we were to prove. 

Theorem 5.3. Let fi denote a bounded simply connected closed region with interior R and 
boundary S. Let u(x,y) be harmonic in R, continuous in 12 and having modulus of continuity 
CO (5) on S. If D is any positive constant, then the modulus of continuity of u{x,y) in 12 satisfies the 
inequality 

co*(6) < o:[D(d/D) 1/3] + (4M/7r) (d/D) 1/3 
for all 8 <D. 

Proof. As in the proof of theorem 5.1, it is sufficient to show that if Cr is a circle 
with center {xi,yi) on S and radius r and if (x2,y2)e^ and {xi—X2y+(yi — y2y<8^<r^, then 

M^2,2/2),C^.nS,12]> l-(4/7r)(5/r)i/l (5.4) 

The theorem then follows if we let r=D{8lDy'^ since 8<D and hence D(8/Dy^^>8. 

Now, if 12 is a Jordan region, (5.4) follows at once from theorem 5.2, since tsin.~^\z\^^^< \z\^^^ 
for 1^1 < 1. We indicate the modifications necessary to include the case where 12 is an arbitrary 
simply connected region. 

Let be a boundary point of 12, let Cr be the circle |^|=r, and let Zq be a point interior to 12 
and to Cr. If no point of 12 lies exterior to Cr, then a function v(z) harmonic in R and continuous 
in 12 equal to unity on S is identically unity in 12 so we have ^(^o) = 1. We proceed to study the 
contrary case. If points of R lie exterior to Cr, such points can be joined to Zq by a Jordan arc 
lying wholly in R, so at least one arc of Cr lies in R. Let the totality of mutually disjoint arcs 
of Cr in R be ^1, A2, .... Denote by Rq the subregion of R interior to Cr containing Zq and 
by Ri the sum of Rq, its reflection (inverse) in Cr, and the arcs Ajc which form part of the boundary 
of Rq] if an arc of Cr is part of the boundary of R, it does not belong to i?i. We modify /?i by 
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adjoining to i?i the interior of each Jordan curve that can be drawn in J?i and to which and Zq 
are exterior. Denote this new region by i?2. Then i?2 contains Zq, is simply connected, and 
has as either an exterior or a boundary point. The subregion Rz of i?2 interior to Or is also 
simply connected, and part of its boundary is an open arc A^ of Or of which every point is 
accessible; we choose A^ as the largest such arc, so that every boundary point of Rz not on Aq is 
either a boundary point of R interior to Gr or is a point of Cr not contained in an arc of Cr 
consisting wholly of accessible boundary points of R^. By a conformal map of Rz onto the 
interior of a circle 7 (in which ^0 necessarily corresponds to an arc of 7) it is clear that H(z,Ao,Rz) 
exists and is unique; this function takes the boundary value unity at every point of Aq and the 
boundary value zero at every boundary point of Rz not on the closure of ^q. Caratheodory^s 
proof of Milloux's theorem is valid for the region i?3, the distance from Zq to the boundary of 
J?3 is not greater than l^ol, so (5.4) is valid for an arbitrary simply connected region, and the 
theorem follows. 

The modulus of continuity of a function F{x,y) on any subset 12^ of L[h,k] is defined by 

0)^(6) =LUB\F(x,,y,) -F{x,,y2)l 

where {xi,yi), {x2,y2) eQ^, and where {xi—X2y+(yi—y2y<^^. 

For d.h. functions in a rectangle we prove the following theorem: 

Theorem 5.4. Let 12 denote the rectangle 0<x<a, 0<y<b, where {a,h)eL[h,k]. Let 12/, 
denote the subset of L[h,k] contained in 12. Ij U{x,y) is d.h. in Rl, the interior of 12x, and has 
modulus of continuity co(6) on Sl, the boundary of Hl, then for b<r={aby^ the modulus of 
continuity of TJ{x,y) in 12/, satisfies the inequality 

i^i^^{b)<o^[2"'{{rby^'+h+k]]+{Ml2)v{\l2){blry" (5.5) 

where M and v{\l2) are defined by (5.2) and (4.6), respectively. 

Proof. As in the proof of theorem 5.1 the maximum value of \U(xi,yi) — U(x2,y2)\ for 
{xi—X2)^+(yi—y2y<^^ is assumed when either (xi,yi) or (2:2,2/2) belongs to Sl. Let us assume 
that {xi,yi)eSL. Now if (0,2,^2) also belongs to Sl the theorem is trivial, since b<r. We there- 
fore assume (0:2,2/2) ei^z,- 

hetp denote a straight line containing {xi,yi) and including one of the sides of the rectangle. 
Let (7 denote the closed interior of a semicircle with center at (xi,yi), with radius bi = 2^^^[ (5r)^/^+ 
h-\-k}, with bounding diameter included in p and containing at least one point of Rl. Since 
b<r we have bi'>2^^^b. Evidently there exists a rectangle T included in C with one side con- 
tained in p and with vertices contained in L[h,k] such that the sides perpendicular to p are 
at least half as long as those parallel to p and such that the latter sides have length at least 
2(5r)i/2. 

By (4.2) and (4.3) we have 

HL[(x,y),CnSLML]>HL[{x,y),TnSL,^L]>HL[(x,y),TnSL,nLf]T]. 

Now, for {x,y)eSLC\ T we have 

HL[{:x,y),T[\vJ\<\=HL[{x,y), TnS^,12^nT]. 

Moreover, since the intersection of y and the boundary of TV[^l is contained in TV\Sl we have 
HL[(x,y), T(]p,T] = for all points of the boundary of ^l^T not contained in Sl^ T. Therefore, 
for all points of the boundary of TH^l we have 

HL[{x,y),TOSL,^LnT]>HL[{x,y),T[]p,T]. 

It follows that for all {x,y)eTr\ilL we have 

HL[(x,y),TnSLML^T]>HL[{x,y),Tnp,T], 
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Now since (x2,2/2)eTni^L7 we have by theorem 4.1 

H,X{x,,y,),mp,T]>l-{ll2)v{lj2){&lry'\ 
Hence, 

Therefore, 

\U{x,,y,)-U{x,,y2)\<o^[2'''{{rby''+h+k]] + {MI^^^^^ 

and the theorem follows. 

In the above proof \ih=^k, one could have obtained a rectangle with the desired properties 
by letting Chave radius 2^'^[{rbY^^-\-h]. Therefore, we have 

Corollary: If h=k, then 

6. Error Estimate for the Rectangle In Terms of the Modulus of Continuity on 

the Boundary 

In tliis section we obtain an error bound for the rectangle of section 3 under the assump- 
tion that the function /(a;,?/), which determines the boundary values, has modulus of con- 
tinuity co(5) on S. The function /a (x,?/) defined in section 3 is also continuous, and we denote 
by co3(6) its modulus of continuity for <x<a, y=0. 

Let Us(x,y) and U3(x,y) be given by (3.3) and (3.5), respectively. We define the function 
^^(x) by the partial sum 

A 

fb{^) = ^dA,nAn sin (riTTx/a), (6.1) 

7J = 1 

where tiie An are the Fourier coefficients for/3(a;,0)=^(;r) as defined by (3.4) and where the 
coefficients dA,n are summation coefficients defined by Jackson [7, p. 9]. Let U5(x,y) and U^{x,y) 
denote respectively the solutions of the Dirichlet problem and its difference analog vanishing 
on S except for y=^. If y=^, the values are determined by/5(ic). We now prove 
Theorem 6.1. IJ (x,7j)eUz. dnd 0<a;<a, e<y<b, then 

\U,(x,y)-u,(x,y)\<SM*CHhlay(e/a)-\ (6.2) 

where 

^=J^+«W^' ^*=Max !,(.)! (6.3) 

and where C and a are determined by (3.10) and (3.12), respectively. 
Proof. It can be verified that 

A 

Ur,{x,y) = '^dA,nAn sin {n'Kxla)Vn{y), (6.4) 

w = l 
A 

U^{Xyy) = YlidA,nAn sin (mrx/a)rm{y), (6.5) 

n = l 

where m and n are related by (3.7), and where Tm(y) is defined by (3.16). Evidently we have 

\U,{x,y)-u,{x,y)\<j:\dA.nAn\'\TUy)-Tn{y)\. (6.6) 

71 = 1 

By lemma 3.3 we have 

^r^(?/)[<(7r?//a) coth {rmrb/a) exp (-mx?//a), {0<y<b), (m>0). (6.7) 
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By lemma 3.2 we have m>an for l<n<A, and 

d 



dm "^-^y^ 



<{'jry/a) coth (arb/a) exp (—amry/a), (6.8) 



Using lemma 3.1 we get, as in the proof of lemma 3.5, 

I r ;;, (^) — r„ (y) I < Cae {hi a) ^ {iryla) v? exp ( — amry/a) . 
Therefore, we have 

T:\^m{y)-T,(y)\<f:^\TUy)-T,{y)\<Cae{h/a)\wyla)J2n' exp [-aniryfa), 

n=l n=l n=l 

which is convergent for y^O. 

Now the following statement can be verified easily: if G{x) is a continuous nonnegative 
function for x>Q, nondecreasing for 0<x<X, and nonincreasing for X<x, then 

^G(n)< rG(t)dt+G(X), 

n=l Jo 

Therefore, since for j^>0 we have 

Max [x^ exp {—vx)] = {3/veyy 
it follows that 

CO /» 00 

y^Ti^ exp (—anTry/a) < t^ exp (—aTryt/a) dt+{3a/aTryey<Q(afairyy-\-{3afaTryey. 
71=1 Jo 

Therefore, we get 

S|^.(^)-^„(2/)|<XC(/l/a)2(y/a)-3<X(7(Va)2(€/a)-^ (6.9) 

n = l 

since y>€. 

For the Fourier coefficients, since \g{t)\<M'^, we have 

|^^| = |(2/a) r git) sin {mrtla)dt\ <{2/a) (I \g{t)\dt<2M'', (6.10) 

Jo Jo J 

Jackson [6] proved that |1— c?^,n|:^37i/A<3. Hence 

\dA,n\<^. (6.11) 

From (6.6), (6.9), (6.10), and (6.11) the theorem now follows. 
Corollary. For (x,y)e^L (^^d 0<x<a, e<y<b we have 

I U^(x,y) -u,{x,y) \<SWC\{hla)'(e/a)-'+2Ko:,{2h), (6.12) 

where K is a positive absolute constant less than ^^ 3 . 

Proof. 

I Us(x,y) —u^{x,y) \ < \ U^{x,y) — U^ix^y) \ + \ U^{x,y) —u^{x,y) \ + \ur,{x,y) — Uz{x,y)\, 
By the maximum and minimum principles for harmonic and d.h. functions, w^e have for (x,y)eQi^ 
I Uz {x,y) — U5 (x,y) 1 , 1^3 (x,y) - u^ (x,y) \ < Max I/5 (x) -g{x)\, 

0<a:<a 

Using a theorem of Jackson [7, p. 7], we obtain 

Max \Ux)-g(x)\<Kc,,{2a/A), 

0<z<a 



10 See Jackson [7, p. 7]. 
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where /^ is a positive absolute constant less than 3. The corollary now follows from the theorem. 
We shall now derive an error bound for the original problem where we no longer assume 
that the boundary values vanish at the corners or along the sides of the given rectangle. 

Theorem 6.2. For all (x,7j)eQL, we have 

|C7(x,2/)-i^(x,^)l<Max|/;2co2[r(My']|, (6.13) 



where 



I=^lQM[\aCaa-'^\a/by^'Kh/ay^' + \,C,a'^\b/ayf'\k/by^']-h 
4K[o:{2h) + c^{2k) + 2M[h/a + klb]}, 

r=(aby^' (6.14) 

^-^ '"'24^/^ coth (a^b/a), 0,="^^^^^^ coth {^^a/b) (6.15) 

a={2/Tra) sinh-^cr), ^^(2(T/Tr) sinh-i(c7-0 (6.16) 

_ 66 27/2 \ 

^" a'TT^'^aVe' ( 

X __:66_ _27/^ ( 

^' pV'^fiVe' ) 



(6.17) 



M=Msixu{x,y) — Mmu(x,y) (6.18) 

(x,y)tS {x,y)eS 

c^t{d) = c^[2'^'{(r8y^'+h + k}]+(M/2)vil/2){dlry^' (6.19) 

and o){8) is the modulus of continuity of u{x,y) on S. The function v{bla) is defined by (4.6) and 
j/(l/2) equals 3.97, approximately. 

Proof. 

Lemma 6.1. If Ui(x,y) is defined by (3.1), then for all (x,y)€Q we have \u(x,y)—Ui(x,y)\<M. 

Proof. Since u(x,y) and Ui(x,y) are harmonic, the maximum of \u(x,y)~Ui(x,y)\ is 
assumed on S. We can assume, without loss of generality that the maximum occurs at a 
point whose ordinate is zero. By (3.1) we have 

Ui(xfi) = (l—x/a)u(Ofl) + (x/a)u(ajO)j 
and 

|u(;r,0)— tii(:r,0)|<|l— :z;/(i|-li6(x,0)— u(0,0)| + |a:/a|.|u(a:,0)-'?^(a,0)|. 
But, 

\u{xfi)— u{Ofl)\,\u(xfi)— u{afi)\< Max u{xfi)— Mm u(xfi)<My 

0<x<a 0<x<a 

and the lemma follows. 

The above result cannot be improved as one can show by considering the case where 
Ui(x,y) is identically zero and u(x,y)>0. 

From lemma 6.1 we conclude that \uz(x,y)\<M. In particular we have \g(x)\<M; hence, 

M*<M. (6.20) 

Lemma 6.2. 

0)3(5) <co(6)+5M/a. (6.21) 

Proof. Evidently we have ooz(d)<o)(8) + ooi(8), where coi(6) denotes the modulus of con- 
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tinuity of Ui(xfl) considered as a function of x. But since 

Ui(x+8,0)-Ui(xfi) = i8/a)[u{afi)—u(Ofl)], 

we have o)i(d)<8M/a, and the lemma follows. 

Lemma 6.3. Ij e<x<a — e,e<y< b — €f{x,y)€^L, then 

\U(x,y)~u(x,y)\<mM[\aCa(h/a)\e/a)-' + \,C,(klb)Xe/by^ 

(6.22) 

Proof. We observe that \U(x,y)—u{x,y)\ = \U2(x,y)—U2(x,y)\. The function [U2(x,y) 
—U2{x,y)] is the sum of [U^(x,y) —U'i{x,y)] and three other terms of similar type. The lemma 
follows from (6.12), from (6.20) and from lemma 6.2. (We note that from the hypothesis 
we hav^e e/a. €/6<K.) 

Lemma 6.3 affords us an error bound for those points of Q^ which are at a distance of not 
less than e from the boundary. In order to obtain a uniform error bound we now consider those 
points which are within e of the boundary. 

Since 4/7r<i/(l/2)/2 it follows from theorem 5.1, (with D=r), and theorem 5.4 that the 
moduli of continuity of U{x,y) and u(x,y) do not exceed the upper bound for cot (8) given by 
theorem 5.4, provided b<{ahy''^. If {x,y)e^L and if {x,y) is within e of a point of S, then there 
exists a point of Sl at a distance not greater than e from {x,y) . We remark that every point of 
Q.L is within {ahy^ of some point of Sl. Since U(x,y) =u(x,y) on Sl we have 

\U{x,y)-u(x,y)\<2c.t(e). 

If / denotes the right member of (6.22), then we have 

\U(x,y)-u{x,y)\<M^x{I'M{e)}, (6.23) 

In order that the error in 12z, should approach zero with the highest power of y=Qiky'^ 
whenever co(5)/5 is bounded as a function of 5, we choose e=r{hk/r^y^'^. Evidently we have 

ela=(j^i\hlaYi\hlayi'\ elb = a-'^'{k/by^'(a/by^'\ 

The theorem now follows from (6.23), lemma 6.3, and theorem 5.4. 

For the special case of the unit square theorem 6.2 and theorem 5.4, corollary, give 
Corollary. If a=b = c=l, then we have 

\TJ{x,y)-u{x,y)\<MdiK[J^]J2], 
where 

J,=^2M\C¥i''+SK[oi{2h) + 2Mhl C={Tr^ll2ae) coth (a7r)~.572, 

J2=^2c^[2''\¥i'-^h)]+Mv{ll2)h'i\ a=(2/7r)sinh-Hl)'-.561, 

_ 66 , 27/2 „ _. 



Substituting numerical values, we obtain 

Ji<66MA2/7+24co(2A)+48MA, J2<2o^[2"\h"''+h)]+4:Mh"\ 

The above expressions for Ji and J2 represent a slight improvement in the formulas previously 
given by the authors [14]. 

7. Other Regions 

In this section we consider regions of the type described in section 2 other than rectangles. 
The case of two overlapping rectangles is studied in detail as an illustration of the method 
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which can be extended to more comphcated regions. We first i)rove a th(H)r(>m vaHd for two 
overlapping regions whether rectangles or not. 

Theorem 7.1. Ltt 12=12'UQ'', where R^{\W is not empty, and where neither region 
includes the other. Let u{x,y) and U{x,y) denote, respectively, the solution oj the Dirichlet problem 
and its finite difference analog Jar 2. Let ^l, ^l, ci't^d Ql denote the points of L[h,k] belonging to 
12, 12', and 12'', respectively. Let Ti^Sin/?!' , 7^1'=^!' ni?L and let 

^, = MaxHr.[(x,y),r,:X], M2=Max ^^[(x,2/),nX]. 

Next, let U' and C7" denote junctions d.h. inQ^ and 12^', respectively, and equal to u on S^ and S^j 
respectively. Then for {x,y)e^L, we have 



where 



\U{x,y)-u{x,y)\<M^x (A,B), 

1— MiM2 {x,y)€Rl 

■M1M2 {x,y)eR'j^ 



^_^ g-(l + M2-MiM2)+g; ij-=Max \U"{x,y)-uix,y)\. 



Proof. Let 



7^' = Max \U{x,y)-u{x,y)\, F^'-Max \U{x,y)-u(x,y)\. 

Clearly, for {x,y)eTL 

\U{x,y)-u{x,y)\<\U{x,y)-W\x,y)\ + \U'\x,y)-u{x,y)\ 

<Mi Max \U{x,y)-W\x,y)\MU"{x,y)-u{x,y)l 

{x,y)tr^ 

Hence, by the maximum and minimum principles for d.h. functions F^<ixiF^^-\-E^\ 
Similarly, F'' <ijl2F'+E\ Therefore, 

j^,^ImE^+E^^ j,,,^ ^2E^^+E^ 

~ 1— jUijLt2 ~" 1— M1M2 

By the maximum and minimum principles, 

Max \Uix,y)-u(x,y)\ <E'+F\ Max \Uix,y)-uix,y)\ <E''+F'\ 

and the theorem follows. 

The proof of this theorem was, of course, motivated by consideration of the difference 
analog of the Schwarz Alternating Process. Indeed, in order to solve the difference equation 
for Ql, one might guess values for U(x,y) on S/CIRl^ and then solve the difference equation 
in Rl obtaining in particular values on Sl'^Rl . Next one solves the difference equation 
for Rl', using the computed values on Sl'^R/, and obtains new values for SlOR/^. This 
process can be repeated and the successive values thus obtained converge to the exact solution 
of the difference equation. Moreover, the rapidity of convergence can be estimated if txi and 
/X2 are known. In fact, one can show that after a complete iteration the maximum error is 
reduced by a factor of (I-/X1JLI2). 

If for 12' and 12" the quantities £'' and £"' are known whenever the modulus of continuity 
and oscillation of u{x,y) are known on /S" and S", then by using theorem 5.1 or theorem 5.3, 
£" and £J" can be computed, provided the modulus of continuity of u{x,y) on S is known. If, 
moreover, jui and/or JLI2 are known, an error bound for the composite region can be obtained. 

We consider now the case of two overlapping rectangles 12i and 122 with sides parallel to the 
coordinate axes. If the intersection of the interiors is not empty, and if neither region is a 
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subset of the other, then the cases A through F shown in figure 1 occur. 



Figure 1. 



The composite regions of cases C and D are of the same type as those of cases A and B, 
respectively. The problem of determining mi and H2 become problems in d.h.m. for the rectangle 
of the types shown in figure 2. 




Upper bounds for d.h.m. are required in points on the closed dotted lines for the open arcs 
on the boundary indicated by heavy lines. 

By (4.3), these upper bounds are not greater than the upper bounds for the problems shown 
in figure 3. 



PROBLEM I 
(FOR TYPE A) 




PROBLEM H 
(FOR TYPE 8.) 



PROBLEM Hr 
(FOR TYPES B2 AND F) 



PROBLEM Ur 
(FOR TYPE E) 



Figure 3. 

For problem I, since the harmonic measure of the open arc is required, ^is an upper bound 
on the dotted line, by symmetry. 

Theorem 7.2. Let U^ be bounded in part by a segment of a line I. Let Rl contain no point 
of I. Let V denote a line perpendicular to I and containing a point of 10 S^. Let S^ denote any 
subset of ICISl contained in one of the open half planes bounded by V . If {x,y) is any point of 
RJ\V, then 

Hd{x,y),S^M<)i' 

Proof. The theorem follows easily by (4.2), (4.3), and problem I. 

As a corollary we have an upper bound of \ for problem II. Also, for problem III if we 
extend the lines whose harmonic measure is required to the right, we get by symmetry as in 
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problem I, ^ as an upper bound for the d.h.m. on the left dotted Hne. Shnilai'Iy ^ is an upper 

bound for the d.h.m. on the right dotted line. 

For problem IV the situation is somewhat more comphcated. We prove 

Theorem 7.3: Let Q. denote the region: x>0, y^O, and let ^l denote the subset of L[h,k] 

contained in Q. Let la, /& denote respectively the intervals 0<Cx^a, y=0 and 0<^y<^b, x — i), where 

the point (a,b) belongs to L[h,k]. We have 



where 



Hd{x,b)JaliIi>,^L]<Max[- 



+w {l'—w)b 



W=(2/Tr) log 



2~'- 



l+exp(— 7r/2) 
1— exp(— 7r/2)* 



22/0 



} 



(7.1) 

(7.2) 



y 




D 


i't 


lb 




B 


Jc(o,b) 


lb 




la 



A E 



Figure 4. 



Here yo is the smallest number not less than a/2a such that yo/k is an integer. The quantities a and c 
are defined by (3.12) and (2.3), respectively. 

Proof. Let /*, /* denote respectively the intervals x>a, y=0, and y^b, x=0. We first 
prove 

Lemma 7.1. If 0<x, 0<y, then 



and 



Hdix,b),I*,Q,] >Hd(x,b),h,Q^], (7.3) 

Hd{a,y),i:,n,] >HA{a,y),IaM- (7-4) 

Proof: By theorem 7.2, we have 

Hi{{x,b),y>h[\x=Q,x>(!i\>Hd{x,h),y<,hr\x=0,x>0]. 

Hence, we certainly have 

i?J(a;,6),{i/>6nx=0}U{-6<2/<Ona;=0},z>0] 

>HjX{x,h),{O<y<b^x=O]\}{y<-h[\x=O],x>0]. 

But by a sign-changing reflection about the hne y=0, we can show that the left member of the 
above expression minus the right member equals i7z,[(a;,6),7j,fii]— Hi[(a;,6),/(,,OJ>0. 

This proves (7.3). Evidently (7.4) can be proved by the same method. 

From lemma 7.1, we conclude that 



77J(x,6),7,U7„fiJ <i/J(^,6),/«,fiJ + (l/2)firj(x,6),0<yna;=0,fid. 



(7.5) 
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But by (7.4) we have 

{ll2)Hd{a,y)fi<y[]x=Q,^r]+Hi{{a,y)Ja,^L]<ll2. 
Therefore, the d.h. function 

W{x,y)=HrX{x,y)Ja,^L] + {ll2)HjX{x,y)fi<y[]x=()M (7.6) 

satisfies the conditions 

Wixfi)^l, {0<x<b); W{0,y)=^^/2, {0<y); W(a,y)<l/2, (0<y), (7.7) 

Let ft* deiiote the region 0<x<af]y^O, and let 12* denote the corresponding subset of 
L[h,k], Evidently, 

W(x,y)<U,{x,y), (7.8) 

where 



(b) L\{xfi) = l, {0<:x<a) ^ (7.9) 



(a) Ui{x,y) is d.h. and bounded in i?* 

(b) L\{xfi) = l, {0<x<a) 

(c) t7i(0,2/) = C7i(a,7/) = l/2, {y>0) 



We now prove 

Lemma 7.2. Let U2(x,y) be d.h. and bounded for (Xjy)eRl and such that 

U2{xfi)=^h (0<x<l); U,{0,y) = UMy)=0, (2/>0). 

We have 

U2(x,c)<Max [w;l-{l-w)c/yol (7.10) 

Proof. Substituting in (3.5) and (3.6) and taking the limit as b becomes infinite we get 

A 

U2(x,y)= S {2/ A) cot {mrl2A) sin (mrx/a) exp {—miry la), 

n odd 

where Ah=a. 

By lemma 3.2, we have m> an. Also, since cot [mrl2A] <2A/mr ior n<A we have 

t^2(x,2/)< S(4M) exp {-Tayn/a)< ±(4.1x71) exp i-irayn/a) = i2M log i"'"!^^ ["''"S ' 
»=i n=i 1— expl^ — Tray/a) 

n odd n odd 



Evidently 


it y^a/2a, we 


have 

Z72(^,^)<^ 












(7.11) 


Now since 0<i^< 1 , 


the function 


















Us{x,y) = l-(l'-w)ylyo 










(7.12) 


is d.h in R* 


and is not less than U2(x,y) for StO y<yo 


and iory= 


=yo. 


Therefore 


for 


y<yo 


we have 






U2{x,y)<l-{1- 


-^)ylyo^ 










(7.13) 
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and the lemma follows. 

Now U2(x,y)=2Ui(x,y)-l. We arlso have by (7.5), (7.6), and (7.8). 

H^[(x,b)Ja[}I„n^] < W(x,b) < U,(x,b), 

and the theorem now follows. 

In applying the theorem to problem IV, we note that we can assume the value zero at the 
point (0,0) without changing the d.h.m. 

Thus, an error bound can be obtained for the case of two overlapping rectangles. The 
above methods can also be used to find error bounds for other regions of the type considered in 
section 2. 
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